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Novel type of atomic-scale spin lattice at a surface and its emergent Hall effect 

M. Hoffmann, 1 ’^ J. Weischenberg, 2 B. Dupe, 1 F. Freimuth, 2 P. Ferriani, 1 Y. Mokrousov, 2 and S. Heinze 1 

institute of Theoretical Physics and Astrophysics, Christian-Albrechts 
University of Kiel, Leibnizstrasse 15, D-2f098 Kiel, Germany 
2 Peter Griinberg Institut and Institute for Advanced Simulation, Forschungszentrum Julich and JARA, Germany 

(Dated: March 13, 2015) 

We predict the occurrence of a novel type of atomic-scale spin lattice in an Fe monolayer on 
the Ir(001) surface. Based on density functional theory calculations we parametrize a spin Hamil¬ 
tonian and solve it numerically using Monte-Carlo simulations. We find the stabilization of a 
three-dimensional spin structure arranged on a (3 x 3) lattice. Despite an almost vanishing to¬ 
tal magnetization we predict the emergence of a large anomalous Hall effect, to which there is a 
significant topological contribution purely due to the real space spin texture at the surface. 


Localized stable spin textures such as skyrmions or 
chiral domain walls have attracted much attention re¬ 
cently due to their unique topological and transport 
properties m and potential applications in spintron- 
ics ms. a key ingredient for their occurrence is the 
Dzyaloshinskii-Moriya (DM) interaction [8j 9i, which 
arises due to spin-orbit interaction in systems with 
broken inversion symmetry, as in the bulk of non- 
centrosymmetric crystals or at surfaces and interfaces. 
Hall effects play an important role in these systems. For 
instance, the spin-orbit torque originating from the spin 
Hall effect drives the motion of chiral domain walls in 
ultrathin films very efficiently and very high speeds have 
been reported 013]. The topological Hall effect, defined 
as the contribution to the Hall resistivity due to chiral 
spin texture, serves as one of the main tools to pinpoint 
the skyrmion phase in the phase diagram of bulk alloys 
such as MnSi or FeGe [ME]. 

The topological Hall effect in complex large-scale mag¬ 
netic structures is normally described assuming the adia¬ 
batic viewpoint of infinitesimally slowly varying spin tex¬ 
ture pQ. For skyrmions, the topological Hall resistivity 
can be factorized into the product of an emergent mag¬ 
netic field, which is the direct consequence of the non¬ 
zero topological charge, and the topological Hall coeffi¬ 
cient _R top , which can be determined from the electronic 
structure of the ferromagnetic crystal pTH fl3] . The va¬ 
lidity of this picture has been demonstrated for large- 
scale skyrmions in bulk Mn^Fei-^Si alloys M- On the 
other side of the length scale, the chirality-driven con¬ 
tribution to the anomalous Hall effect (AHE) has been 
predicted and observed in bulk strongly-frustrated cor¬ 
related oxides and bulk antiferromagnets, which exhibit 
non-collinear magnetic order on the scale of 1 nm [T5ll22j . 

Chiral domain walls and skyrmions with an extent of 
down to 1 nm can also occur at transition-metal inter¬ 
faces and surfaces j23b26] . It has been demonstrated that 
higher-order exchange interactions can play a crucial role 
in such systems. They enforce the atomic-scale skyrmion 
lattice observed for an Fe monolayer on the Ir(lll) sur¬ 
face [25] and a conical spin spiral phase for a Mn double 
layer on W(110) [27] . However, to date, very little is 


known both experimentally and theoretically about the 
Hall effects in such complex nanometer-scale spin tex¬ 
tures at surfaces and interfaces. 

Here, we present a model system for a novel type of 
atomic-scale spin lattice at a transition-metal surface and 
study its transport properties. In contrast to systems 
explored so far for skyrmion spin textures the local ex¬ 
change interaction is antiferromagnetic in this system. 
Based on density-functional theory (DFT) and Monte- 
Carlo simulations, we find a complex three-dimensional 
spin structure on a (3x3) lattice for an Fe monolayer on 
Ir(001), which has been shown to grow pseudomorphi- 
cally [28]. This spin texture with angles close to 120° 
between adjacent spins arises due to the coupling of spin 
spirals stabilized by exchange and DM interaction via the 
four-spin interaction. Using DFT we demonstrate that, 
despite very small total magnetization, this non-trivial 
(3x3) spin lattice gives rise to the AHE, essentially larger 
than that observed in non-collinear bulk compounds H3- 
m. We further show the existence in this system of a 
sizeable surface “topological” contribution to the AHE 
and orbital magnetization, which originate purely in the 
non-trivial real-space spin distribution and do not rely 
on the presence of spin-orbit interaction. 

Nanoscale spin textures at transition-metal inter¬ 
faces Eaunu can be treated employing a Hamiltonian 
on the discrete atomic lattice 
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which describes the magnetic interactions between the 
magnetic moments M/ of atoms at sites R/. For the Fe 
monolayer on Ir(001), we used DFT to obtain the pa¬ 
rameters for the exchange interaction (Jr/), the DM in¬ 
teraction (D/j), the four-spin interaction [K^ki) and the 
biquadratic exchange ( Bij ) as well as a uniaxial magne- 
tocrystalline anisotropy (K±). We applied the projector 
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FIG. 1: Energy dispersion of homogeneous flat spin spi¬ 
rals for Fe/Ir(001). The energies E{ q) (filled circles) are 
calculated via DFT along the high symmetry lines of the 
two-dimensional Brillouin zone and given with respect to the 
c(2 x 2) antiferromagnetic state. The solid lines are fits to 
the Heisenberg model with up to sixth nearest-neighbors [36]. 
The inset shows the energy dispersion close to the M point for 
left- and right-rotating spirals including SOC, i.e. the effect 
of the DM interaction. 


augmented wave (PAW) method [32] as implemented in 
the VASP code [53J [34]. Computational details are given 
in the supplementary material. 

To determine the exchange constants we have 
considered flat spin spirals in which the magnetic mo¬ 
ments are confined in a plane with a constant angle 
between moments at adjacent lattice sites propagating 
along high symmetry directions of the surface. Such a 
spin spiral can be characterized by a wave vector q from 
the two-dimensional Brillouin zone (BZ) and the mag¬ 
netic moment of an atom at site R^, given by = 
M(sin (qRi), cos (qRi), 0) with the size of the magnetic 
moment M. 

The calculated energy dispersion E( q) of spin spirals 
for Fe/Ir(001) is displayed in Fig. [T] At the high symme¬ 
try points we obtain collinear spin structures: the fer¬ 
romagnetic state at T, the c(2 x 2) antiferromagnetic 
state at M, and the p( 2 x 1) antiferromagnetic state at 
X. Clearly, the c(2 x 2) antiferromagnetic state is low¬ 
est in energy among the considered collinear states in 
agreement with previous DFT studies [35] • The energy 
dispersion is very flat in the vicinity of the M-point due 
to the frustration of exchange interactions. A fit to the 
Heisenberg model, i.e. the first term in Eq. 0- with Jij’a 
up to sixth nearest neighbors [36] leads to an excellent 
description as shown by the solid line in Fig. 

Note, that the energy dispersion of Fe/Ir(001) is almost 
inverted with respect to Fe/Ir(lll) where the energy dis¬ 
persion is flat around the f-point, i.e. the ferromagnetic 
state O EH- Therefore, we can also expect complex 
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FIG. 2: Spin structures of lowest energy obtained via Monte- 
Carlo simulations as a function of the strength of the four-spin 
interaction. The biquadratic interaction is changed according 
to the condition 2iF4 sp i n + B — 0.7 meV from the DFT cal¬ 
culations. The sketches only display a small section of the 
actually simulated spin lattice. 


three-dimensional spin structures to occur here but of 
different type due to the nearest-neighbor antiferromag¬ 
netic exchange. 

By taking spin-orbit coupling (SOC) into account, we 
can determine the magnetocrystalline anisotropy energy 
(MAE) defined as the energy difference between config¬ 
urations with different orientation of the magnetization. 
For the collinear state of lowest energy, i.e. the c(2 x 2) 
antiferromagnetic state, we found an easy out-of-plane 
axis with a MAE of K± = —0.25 meV. 

At a surface SOC also induces the DM interaction [39, 
00]. In order to determine its strength, we have calcu¬ 
lated the total energy of a 120° spin spiral along the 
TM-direction in a (3 x 1) super cell including SOC both 
with a left-handed and a right-handed rotational sense. 
We find that spin spirals with a right-handed rotational 
sense are lower by 7.3 meV/Fe atom. This energy differ¬ 
ence allows to calculate the value of the DM interaction 
within the nearest-neighbor approximation which results 
in a value of Di = 1.5 meV. Including the DM inter¬ 
action into the energy dispersion of spin spirals leads to 
an energy minimum at an angle of about 138° between 
adjacent spins as shown in the inset of Fig. [l] 

From the energy dispersion of spin spirals, only the 
Heisenberg-type exchange can be obtained. The impact 
of higher-order spin interactions can be determined by 
considering superposition states of two spin spirals. If 
only Heisenberg-type exchange played a role all of these 
spin states would be degenerate in energy. However, our 
DFT calculations show considerable energy differences 
on the order of a few meV/Fe-atom (see supplementary 
material). From these calculations, we determine that 
the nearest-neighbor four-spin, iC 4 sp i n , and biquadratic, 
H, interaction fulfill the condition 2 K 4 sp [ n +B = 0.7 meV. 

The energy functional Eq. 0 with the parameters 
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FIG. 3: (a) Side view and (b) top view of the proposed 

atomic-scale (3x3) spin lattice. Two unit cells are shown. 


from DFT can be minimized using Monte-Carlo simu¬ 
lations based on the Metropolis algorithm. We have cho¬ 
sen a spin lattice of (66 x 66 ) spins and used periodic 
boundary conditions. We have checked the impact of the 
lattice size and of using open boundary conditions and 
found no effect on the obtained ground state. In order 
to explore the impact of the higher-order spin interac¬ 
tions, which are not univocally determined by our DFT 
calculations as discussed above, we have chosen different 
values of B and K 4 sp [ n that are in accordance with the 
condition given above. We changed the value of the four- 
spin interaction in steps of 0.1 meV and the biquadratic 
interaction and J 3 were modified accordingly [36]. 

We found three different types of ground states de¬ 
pending on the value of FQgpm as shown in Fig. [ 2 } A 
large biquadratic interaction results in a so-called up- 
up-down-down ( uudd ) state since a collinear alignment 
of neighboring spins is preferred. However, if the bi¬ 
quadratic interaction is reduced we find an atomic-scale 
non-collinear (3x3) spin lattice that is stabilized by the 
four-spin term. For a value of K 4 sp [ n > —0.4 meV the 
four-spin term cannot couple the spin spirals and we ob¬ 
tain a spin spiral ground state with an angle of about 
140° between adjacent spins. 

As shown in Fig. [ 2 J the novel (3x3) spin lattice can 
occur for a large range of values of the four-spin interac¬ 
tion. Its spin structure is shown in Fig. [3j The spins at 
the corners of the unit cell point upwards perpendicular 
to the surface while the spins along the sides rotate with 
an angle of ~ 123° from the surface normal. The four 
spins in the center of the cell point towards the corners 
and with an angle of ~ 22 ° out of the film plane. 

The occurrence of this three-dimensional spin struc¬ 
ture can be understood from the interplay of the differ¬ 
ent interactions. The combination of exchange and DM 
interaction leads to a spin spiral with an angle of approx¬ 
imately 120 ° between adjacent spins and thus a periodic¬ 


ity of 3 atoms (cf. Fig. [T]) . For Fe biatomic chains on the 
(5 x 1 ) reconstructed Ir( 001 ) surface, such a spin spiral 
state has been experimentally observed [24]. In the Fe 
monolayer on Ir( 001 ), the four-spin interaction can cou¬ 
ple these spin spirals into a square lattice. Note, that 
there is an opposite rotational sense of the spin rotation 
along the side and the diagonal of the unit cell. This 
results from the antiferromagnetic exchange coupling be¬ 
tween nearest neighbors which is stronger than the DM 
interaction that would prefer a unique rotational sense 
along both directions. Due to this peculiar competition 
of DMI and Heisenberg exchange, the spin lattice is ex¬ 
tremely stable in an external magnetic field and cannot 
be destroyed up to 80 T as found in our MC simulations. 
The transition temperature to the paramagnetic state is 
obtained at approximately 60 K. 

To investigate whether the (3x3) spin texture results 
in non-trivial transport properties, we compute from first 
principles m the intrinsic Berry curvature contribution 
to the ^-component of the anomalous Hall conductivity 
(AHC) in the system cr^g = J BZ ^^(k) dk [42] . 
where 


n xy (k)= E 2Im 

ti<Ef m^n 


(V^nk \Vx | V4nk) (V^mk \Vy | V^nk) 
(^rak ^nk)^ 


( 2 ) 

is the Berry curvature of occupied states with ^ n k as 
the Bloch states with corresponding energies £ n k, and Vi 
is the i’th Cartesian component of the velocity opera¬ 
tor. The results of our calculations for cr^ 3 , presented 
in Fig.[4]as a function of the substrate thickness, indicate 
a sizeable AHE in the (3x3) spin lattice state with the 
magnitude similar to that of bulk transition-metal ferro- 
magnets (42H44] . The large variation of the AHC with 
thickness, apparent from Fig. [4] is typical for such effects 
as the AHE, spin Hall effect or the spin-orbit torque in 
the limit of ultrathin films [43,. 45, [46]. 

In the context of thin magnetic layers of transition- 
metals on paramagnetic substrates, the emergence of the 
large cr ^ 3 appears rather surprising, since the total mag¬ 
netization of the system in the (3x3) state is almost van¬ 
ishing. By artificially rotating the spin moments on the 
Fe atoms slightly away from their equilibrium directions 
we acquire a complete suppression of the magnetization 
and observe that the values of < 7^*3 stay very close to 
those with small uncompensated magnetization. This 
clearly distinguishes our case from the case of the AHE in 
collinear magnets, which relies on non-vanishing macro¬ 
scopic magnetization and presence of SOC [42J. 

Another remarkable observation is that a large con¬ 
tribution to cr ^3 is provided even without taking the 
SOC into account, as apparent from Fig. [4] where the 
values of the intrinsic AHC, computed with the SOC 
explicitly switched off in our calculations, are presented 
in comparison with cr^g. Since the AHE vanishes for 
any collinear magnetic state of our system without SOC, 
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FIG. 4: Top: BZ distribution of the Berry curvature without 
(left) and with (right) SOC for Fe monolayer in (3 x 3) state 
with one layer of Ir substrate, superimposed with the real- 
space distribution of the spins (blue arrows). Bottom: calcu¬ 
lated values of (T^x 3 as function of the Ir substrate thickness. 

it allows us to define the contribution to < 7^*3 without 
SOC as the “topological” contribution to the AHC, cr™ 3 , 
which stems purely from the spin texture in real space, 
and which does not rely on the presence of SOC. The 
particular symmetry of our system which results in non¬ 
vanishing cr™ 3 , also gives rise to a finite local scalar spin 
chirality • (Mj x M&), non-vanishing when integrated 
over the unit cell m- To distinguish our case from the 
case of large two-dimensional skyrmions and bulk frus¬ 
trated oxides, for which topological contribution to the 
Hall effect in some cases can be described neglecting the 
spin-orbit effects mmmm, for our class of systems 
we call the corresponding anomalous Hall effect without 
SOC the surface topological Hall effect. Our calculations 
suggest the existence of surface THE in transition-metal 
multilayers. 

Ultimately, the large values of cr ™ 3 are due to a direct 
effect of the non-trivial real-space distribution of spin on 
reciprocal-space distribution of the AHC, given by the 
Berry curvature m- To convince ourselves in this ex¬ 
plicitly, we plot in Fig. [4] the Brillouin zone distribution 
of the Berry curvature computed with and without SOC 
for the system of an Fe layer in the (3 x 3) spin state on 
one layer of the Ir substrate. As apparent from the case 
without SOC, there is a very close correlation of the Berry 
curvature distribution with the spin-distribution in real- 
space, while the effect of SOC is to provide an additional 
fine structure to this distribution stemming from SOC- 
induced band splittings in the vicinity of the Fermi level. 
Thus the surface THE is more complex than the THE in 
large-scale skyrmions for which the topological contribu¬ 
tion to the THE — the emergent field — can be separated 
from the electronic effects in a collinear host encoded in 
R top QHl 03] . The surface THE arises from a close inter¬ 


twining of the real and reciprocal space topology, which 
together play a role of a single multi-dimensional topo¬ 
logical object with non-trivial transport properties. 

The microscopic origin of the competition between 
non-collinearity and spin-orbit interaction for the AHE in 
such non-trivial surface spin textures as considered here, 
presents an exciting direction to study both theoretically 
and experimentally [48J09]. In particular, we conjecture 
that the surface THE is commonly an important part of 
the AHE exhibited by complex spin structures at sur¬ 
faces, such as nanoskyrmions [25]. One of its prominent 
manifestations would be the contribution to the orbital 
magnetization at the surface which is not originated in 
spin-orbit coupling mug. The orbital magnetization 
and the Hall effect have the same symmetry and indeed, 
our calculations reveal the formation of non-vanishing 
local atomic orbital moments at the surface of our sys¬ 
tem without spin-orbit. Without SOC, the magnitude of 
the maximal local orbital moment among the Fe atoms 
ranges from — 0.13/i# to 0.03 /ib depending on the sub¬ 
strate thickness, which is similar in magnitude to that ob¬ 
tained with spin-orbit interaction included. Such “topo¬ 
logical” orbital magnetization could be readily addressed 
experimentally by surface techniques. 
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